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on  Pur^  and  Applied  Mathematics, 


Further  Work 

1)   Propagation  of  Shock  Waves  Through  Non-uniform  Regions, 
(H.  B,  Keller.  D,  A.  Levine,  D,  Loveland) 
The  numerical  procedure  developed  in  the  paper  "Motion  of 
a  bore  over  a  sloping  beach"  has  been  extended  for  application 
to  the  problem  of  shocks  moving  through  non-uniform  ducts.   The 
finite  difference  equations  employ  a  fixed  Eulerian  net  in  the 
(x,t)  plane  and  an  iterative  method  for  fitting  in  the  shock  at 
each  time  increment.   This  method  has  many  advantages  over 
previously  employed  characteristic  net  calculations.   In  addition, 
the  stability  and  convergence  of  the  method  have  been  demonstrated 
and  sufficient  conditions  for  the  convergence  of  the  shock-fitting 
iterations  have  been  determined, 

A  general  code  which  employs  this  numerical  procedure  has 
been  written  for  the  I,B,M,  70i|,   The  channel  area  in  this  program 
can  be  of  the  general  fo-rm: 

A(x)  =1,   X  <  Oj   A(x)  =  arbitrary  function  of  x,  x  >  0  , 

Initially  (at  t  =  0)  a  shock  can  be  at  any  position  x  =  £,{0)   and 
is  assiomed  propagating  in  the  positive  x-direction  with  a  uniform 
gas  ahead  of  the  shock  (i.e.  in  x  >  <?(0)),   The  gas  is  assiomed 
polytropic  with  arbitrary  y  (>  1), 

Test  calculations  are  being  performed  at  the  present  time. 
When  the  code  is  completely  debugged  it  is  planned  that  a  series 
of  calculations  with  varying  channel  shapes  will  be  run.   If  the 
results  show  significant  improvement  over  previous  calculations 
or  if  novel  shapes  are  included  a  paper  v;ill  be  written  describing 


,  ■'■'■'.> 


r-:    ^];»0-". 


nc 


{ ;■ . .T'lo-'Ci;  »•'!  .oiixvi^J  «iv  ,u  .--.r.r 


W-.'iS 


\f 


this  material  in  detail.   Due  credit  will  be  given  to  this 
contract  in  any  subsequent  publications  related  to  this  work, 

2)  Structure  of  Magnetohydro dynamic  Shock  Waves,  (G,  B,  Whitham) 

The  discussion  of  this  topic  contained  in  "Some  comments 
on  wave  propagation  and  shock  wave  structure"  (G,  B,  Whitham, 
Comra,  Pure  Appl,  Math,,  Feb,  1959)  has  been  extended  to  include 
a  normal  magnetic  field  and  oblique  shocks.   One  interesting 
result  is  the  non-imiqueness  of  the  shock  profile  in  certain 
cases.   This  result  was  noted  previously  by  G,  S,  S,  Ludford  in 
"The  structure  of  a  hydromagnetic  shock  in  steady  plane  motion" 
(J.P.M, ,  1959).  Ludford  claims  that  all  but  one  of  the  profiles 
can  be  eliminated  as  being  vinrealistic.  However,  we  claim  that 
his  arguments  are  wrong  and  that  the  lack  of  uniqueness  is 
genuine.   Plausible  arguments  are  provided  to  explain  why  the 
profile  is  not  unique  in  these  cases. 

Some  further  results  have  been  obtained  on  this  problem. 
However,  it  has  recently  been  discovered  that  P,  Germain  has 
considered  such  problems  and  previously  obtained  some  of  our 
results.   Thus  it  is  not  clear  that  a  paper  should  be  written  on 
this  work  but  due  credit  will  be  given  to  this  contract  in  the 
event  that  such  a  paper  is  published, 

3)  Thjpee-shock  Intersection,   (Akira  Sakurai,  J,  Hesse,  H,  B, 
Keller) 

Dr,  Sakurai  has  proposed  new  ideas  to  resolve  the  discrepancy 
between  the  simple  three-shock  theory  of  Mach  reflection  and 
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experiments.   He  Introduces  viscous  effects  in  the  sector  which 
contains  the  vortex  sheet,  on  the  grounds  that  the  discontinuities 
are  weak  in  the  cases  where  the  theory  does  not  check, 

A  rather  complete  numerical  solution  of  Dr,  Salmral's  theory 
has  been  obtained.   These  results  indicate  better  agreement  with 
the  experiments  than  does  the  simple  three-shock  theory.   This 
work  has  been  described  in  an  vinpublished  paper.   There  is  still 
some  doubt  as  to  the  validity  of  the  theory  and  until  these 
doubts  can  be  eliminated  no  report  will  be  issued.   If  such  a 
paper  or  report  is  published  this  contract  will  receive  due 
credit, 

i|)   Boundary-layer  Transition  to  Turbulence,   (D,  F,  De  Santo, 
H.  B,  Keller) 

This  work  is  a  rather  complete  revision  of  some  studies 
initiated  by  H.  P,  Ludloff  under  A,P,0,S,R,  -  A,R.D,C,  Contract 
AP  18 (603) -25,   The  problem  is  to  derive  from  the  time  dependent, 
viscous  equations  of  motion  some  of  the  phenomena  that  occur 
during  the  transition  from  laminar  to  turbulent  flow  along  a 
flat  plate.  It  is  assumed  that  well  within  the  boundary  layer 
a  small  harmonic  disturbance  is  present  and  the  deviations  from 
the  Blasius  flow  are  sought.  However,  no  linearizations  are 
employed.  In  terms  of  the  vorticity,  U)  ,  and  stream  function,  cj), 
the  Navier-Stokes  equations  reduce  to  a  non-linear  parabolic 
partial  differential  equation  coupled  with  a  Poisson  equation. 
The  domain  of  interest  is  chosen  as  some  sufficiently  large 
rectangle  downstream  from  the  disturbance. 
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In  this  rectangle  the  equations  of  motion  are  approximated 
by  appropriate  difference  equations  and  numerical  solutions  are 
sought.   The  difference  equations  have  been  shown  to  be  con- 
vergent and  stable,  under  a  mild  restriction  on  the  net  spacing, 
(It  is  believed  that  the  method  is  \mconditionally  stable,  but 
this  has  not  yet  been  demonstrated,)   The  difference  form  of  the 
Poisson  equation  is  solved  by  a  rapidly  converging  iteration 
while  the  semi-implicit  parabolic  difference  equations  are  solved 
explicitly. 

Some  test  calculations  have  been  performed  for  small 
amplitude  (almost  sonic)  disturbances.  However,  a  high  frequency 
wave  develops  and  grows,  and  its  cause  has  been  traced  to  the 
treatment  of  the  boundary  conditions  on  the  upstream  boundary  of 
the  rectangle.   At  the  present  this  difficulty  is  being  investi- 
gated and  some  methods  to  eliminate  it  have  been  proposed. 
Further  calculations  are  planned  and  this  contract  will  be 
credited  in  any  subsequent  publications  resulting  from  this  work, 

5)   Solar  Spicules,   (H,  B.  Keller,  R,  Lilst,  A,  Lapidus) 

This  problem  is  concerned  with  the  formation  and  focusing 
of  magnetohydro dynamic  shocks  in  a  solar  chromosphere.   One  of 
the  basic  objectives  is  to  devise  and  test  numerical  procedures 
for  treating  such  problems  in  two  space  dimensions.  In  the 
particular  problem  under  study  a  perfectly  conducting  gas 
(Lundquist  model)  initially  in  equilibrium  under  gravity,  with  a 
vertical  magnetic  field  present,  is  subjected  to  a  transient 
localized  (pressure)  pulse  near  the  bottom  surface  (assviraed  rigid). 
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The  nvimerical  procediore  used  Is  a  modification  of  Lax's 
conservation  form  for  the  difference  equations.  Artificial 
boTondarles  must  be  Introduced,  due  to  the  limitations  of  existing 
computing  equipment,  and  have  the  effect  of  causing  some  reflec- 
tions of  any  waves  incident  on  them.   The  treatment  at  these 
boundaries  has  been  modified  several  times  in  order  to  reduce 
the  effect  of  such  reflections,   l\/hen  the  magnetic  field  is 
absent  (i,e,  ordinary  hydrodynamics)  the  procedvire  works  quite 
well  showing  the  formation  of  well  defined  shocks  which  propagate 
as  expected.   If  gravity  is  absent  the  results  are  even  better 
(since  uniform  orders  of  magnitude  are  then  maintained  in  the 
ambient  data)  and  comparisons  can  be  made  with  accurate  (but 
approximate)  analytical  solutions. 

However,  when  both  magnetic  field  and  gravitational  effects 
are  present  some  difficulties  remain.   Shocks  can  be  seen  to  form 
and  are  even  focused  along  the  field  but  spurious  disturbances 
still  arise  from  the  top  boundary.   These  signals  grow  as  they 
propagate  downward;  this  is  to  be  expected  for  the  magneto- 
hydrodynaralc  "fast  waves".   Whether  these  spurious  effects  can 
be  tolerated  or  not  depends  upon  their  relative  amplitude  and  how 
far  they  propagate  before  the  relevant  phenomena  have  occurred. 
Methods  for  reducing  these  disturbances  are  being  considered  and 
further  computations  are  in  progress, 

6)   Secondary  Shocks  Vfithln  Spherical  and  Cylindrical  Blasts, 
(M,  P,  Friedman) 

Investigations  into  the  behavior  of  the  gas  flow  behind 
spherical  or  cylindrical  blasts  have  shown  the  presence  of  a 
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secondary  shock.   This  secondary  shock,  at  first  weak,  is 
initially  carried  outward  with  the  expanding  gases;  as  it  develops 
it  strengthens  and  then  bends  back  toward  the  origin  arriving 
there  with  great  intensity.   Descriptions  of  the  flow  field  due 
to  a  spherical  or  cylindrical  detonation  can  be  obtained  by 
nijitierically  integrating  the  complete  equations  of  motion  (H,  L, 
Erode,  Phys,  of  PI,,  April  1959),   However,  using  some  recently 
developed  techniques  in  shock  dynamics  (R,  F,  Chisnell,  J,  Fl, 
Mech, ,  1957;  G,  B,  Whitham,  Comm.  Pure  and  Appl,  Math,,  1952)  it 
is  possible  to  obtain  explicit  formulae  describing  the  main  blast 
wave  and  the  formation  and  subsequent  motion  of  the  secondary 
shock.   The  main  blast  wave  is  described  by  using  Chisnell 's 
shock-area  rule  while  the  secondary  shock  is  obtained  initially 
by  the  VJhitham  weak  shock  approximation.  When  the  secondary 
shock  is  strengthened  beyond  the  region  of  validity  of  the  weak 
shock  approximation  a  generalized  shock-area  rule  is  used,  the 
generalization  being  in  adapting  it  to  the  case  of  a  non-uniform 
flow  ahead  of  the  shock, 

4  4* 

These  formulae  are  now  being  programmed  for  the  I,B,M,  ']o\\. 
computer  and  a  test  problem,  which  has  been  solved  by  numerical 
integration  of  the  flow  equations  (Erode)  will  be  used  to  check 
out  the  present  theory.  Any  subsequent  publications  resulting 
from  this  work  will  be  credited  to  this  contract. 
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Chicago  1,    Illinois 


(1) 


University  of  Maryland 
Institute  for  Fluid 
Dynamics  and  Applied  Mathematics 
College  Park,  Maryland 

University  of  Maryland      (1) 
ATTN:   Engineering  Library 
College  Park,  Maryland 

Massachusetts  Institute  of   (1) 

Teclinology 
Naval  Supersonic  Laboratory 
Cambridge  39,  Massachusetts 

Massachusetts  Institute  of   (1) 

Technology 
ATTN:   Library 
(Please  route  to  Mech,  and 

Aero.  Engr.  and  Mechanics) 
Cambridge  39,  Miassachusetts 

Massachusetts  Institute  of   (1) 

Technology 
Fluid  Dynamics  Research  Group 
)   ATTN:   Dr.  Leon  Trilling 
Cambridge  39,  ^lassachusetts 

University  of  Michigan      (l) 
Department  of  Aeronautical 

Engineering 
ATTN:   Library 
East  Engineering  Building 
Ann  Arbor,  Michigan 
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Midwest   Research   Institute         (1) 

iTTN:      Library 

|25  Volker  Boulevard 

Kansas   City  10,    Missouri 

Jniversity   of  Minnesota  (1) 

Institute   of  Technology 
iTTN:      Engineering  Library 
Minneapolis,    Minnesota 

^osemount  Aeronautical  Laboratories 
Jniversity   of  Minnesota 
Department   of  Aeronautical 

Engineering 
iTTN:      Library 
ilnneapolis,    Minnesota  (1) 

^orth  Carolina   State   College      (1) 
Division   of  iingineerlng  Research 
iTTN:      Technical   Library 
Raleigh,   North  Carolina 

Mo  State   University  (1) 

Department   of  Aeronautical 

Engineering 
iTTN:      Library 
lolumbus,    Ohio 

i'olytechnic    Institute    of  Brooklyn 

Department   of  Aeronautical  Engin- 

iering  and  Applied  Mechanics 

lTTN:      Library 

i33  Jay   Street 

Srooklyn  1 ,    New  York  ( 1 ) 

lerodynamics   Laboratory  (1) 

polytechnic    Institute    of  Brooklyn 
)27  Atlantic  Avenue 
''reeport.    New  York 

Pennsylvania   State   University      (1) 
)epartment   of  Aeronautical 

Engineering 
tTTN;      Library 
Jniversity   Park,    Pennsylvania 

'he  James   Forres tal  Research  Center 
^inceton  University 
^TTN:      Library 

Please   route    to    Prof,    S,    Bogdonoff) 
¥inceton.   New  Jersey  (1) 
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Princeton  University       (1) 
Department  of  Aeronautical 

Engineering 
ATTN:   Library 
Princeton,  New  Jersey 

Rensselaer  Polytechnic  Institute 
Department  of 

Aeronautical  Engineering 
ATTN:   Library 
Troy,  New  York  (l) 

University  of  Southern  California 

Engineering  Center 

ATTN:   Library 

35l8  University  Avenue 

Los  Angeles  1,   California   (1) 

Stanford  Research  Institute  (1) 
Documents  Center 
ATTN:   Acquisitions 
Menlo  Park,  California 

Stanford  University        (1) 
Department  of 

Aeronautical  Engineering 
ATTN:   Library 
Stanford,  California 

Defense   Research  Laboratory 

University  of  Texas 

Post  Office  Box  8029 

Austin  12,  Texas  (1) 

University  of  Virginia     (1) 
Ordnance  Research  Laboratory 
Charlottsville,  Virginia 

University  of  Washington   (1) 
Department  of 

Aeronautical  Engineering 
ATTN:   Library 
Seattle,  V/ashington 

New  York  University        (1) 
Institute  of  Mathematical  Sciences 
ATTN:   Library 
New  York  3*  New  York 


r  I 


ecfc 


,  ■  >'-v,-o.:.!;'r. 


;0.t 


lOiV---^' 


(X) 


.  i  t  v.- 


INDUSTRIAL   0:.GAr'TZjiTIONS 


Allied  Research  Associates      (1) 

ATTN:      Library 

(Please   route   to  Dr.    T,   R.    Goodman) 

14.3  Leon  Street 

Boston  5j  Massachusetts 

AVCO  Manufacturing  Company   (l) 
Research  Laboratories 
ATTN:   Chief,  Technical  Library 
2385  Revere  Beach  Parlcway 
Everett  [|.9,  Massachusetts 

AVCO  Manufacturing  Company   (1) 
Research  and  Advanced  Development 

Division 
ATTN;   Research  Library, 

Mrs.  H,  M,  Page 
201  Lowell  Street 
Wilmington,  Massachusetts 

Bell  Aircraft  Corporation    (1) 
ATTN:   Library 
Post  Office  Box  1 
Buffalo  5,  NeiNf  York 


Boeing  Airplane  Company 
ATTN:   Library 
Post  Office  Box  3107 
Seattle  II4.,  Washington 


(1) 


Chance -Vought  Aircraft,    Inc.       (1) 
ATTN:      Library 
Dallas,   Texas 


CONViilR 

Port  Worth  Division 
ATTN:      Library 
Fort  Worth  1,    Texas 

CONVAIR 

ATTN :      Library 

Post    Office   Box  1011 

Pomona,   California 


(1) 


(1) 


CONVAIR  (1) 

Scientific  Research  Laboratory 

ATTN:   Library 

Post  Office  Box  950 

San  Diego  12,  California 


GONVaIR  (1) 

Scientific  Research  Laboratory 
ATTN :   Library 
(Please  route  to  Chief, 

Applied  Research) 
Post  Office  Box  9^0 
San  Diego  12,  California 

Cornell  aeronautical       (1) 

Laboratories,  Inc. 
ATTN:   Library 
11)4.55  Genesee  Street 
Buffalo  21,  New  York 

Douglas  Aircraft  Company,  Inc. 

ATTN :   Library 

827  Lapham  Street 

El  Segundo,  California     (1) 

Douglas  Aircraft  Company,  Inc. 
ATTN:   Library 
3000  Ocean  Park  Boulevard 
Santa  Monica,  California   (1) 


Pairchild  Engine  and 
Aircraft  Company 
Guided  Missiles  Division 
ATTN:   Library 
Wyandanch,  L.I.,  New  York 

Flight  Sciences  Laboratory 
ATTN:   Library 
1965  Sheridan  Avenue 
Buffalo  23,  New  York 


(1) 


(1) 


(1) 


General  Applied 
Science    Laboratories,    Inc. 
Meadowbrook  National   Bank  Building 
60  Hempstead  Avenue 
Hempstead,    New  York 

General  Electric   Company        (1) 
Aircraft   Gas   Turbine   Division 
ATTN;      Library 
Cincinnati  l5,    Ohio 

General  Electric   Company        (1) 

Aeroscience   Laboratory   -  MSVD 

ATTN:      Library 

(Please   route   to  Dr.    H,   Lew) 

3750   "D"    Street 

Philadelphia  2I|.,  Pennsylvania 


(X; 


leneral  Electric   Company 
lesearch  Laboratory 
>ost   Office   Box  1088 
ichenectady  S,    New  York 


(1) 


eneral  Electric   Company  (1) 

pecial  Defense    Products  Division 
198  Chestnut   Street 
'hiladelphia  h,,    Pennsylvania 

ruiranari  Aircraft  (1) 

ngineering  Corporation 

TTN:      Library 

ethpage ,  L.i,,  New  York 

ughes  Aircraft  Company     (1) 
esearch  and  Development 

Laboratories 
TTN:   Library 
ulver  City,  California 

ockhead  Aircraft  Corporation   (1) 

TTN:      Library 

ost   Office   Box  55l 

urbank,    California 

ockheed  Aircraft  (1) 

issile   Systems  Division 

TTN:      Library 

alo  Alto,    California 

he  Martin  Company  (1) 

TTN:      Library 
altimore   J,,    Mary_and 

cDonnell  Aircraft   Corooration 

TTN:      Library 

ost   Office   Box  5l6 

t,   Louis  66,   IHssourl  (1) 

orth  American  Aviation,    Inc.      (1) 

issile   Division 

FTN :      Library 

22llj.  Lakewood  Boulevard 

Dwney,   California 

srthrop  Aircraft,    Inc.  (1) 

PTN:      Library 
iwthorne,   California 

Lasmadyne   Corporation  (1) 

rTN :  Library 

339  South  I%in  Street 

mta  Ana,  California 


The  Ramo-V/ooldridge  Corporation 

ATTN:   Chief  Librarian 

5730  Arbor  Vitae 

Los  Angeles  l^S ,   California   (1) 


Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California 


(1) 


Republic  Aviation  Corporation 

ATTN:      Library 

Farraliigdalo ,   LoI,,.-New  York      (1) 


RIAS,  Inc.  (1) 

ATTN:   Library 
7212  Bellona  Avenue 
Baltimore  12,  Maryland 

United  Aircraft  Corporation   (1) 

Research  Depart;ment 

ATTN:   Library 

[{.00  Main  Street 

East  Hartford  8,  Connecticut 

VITRO  Laboratories  (i) 

V/est   Orange   Laboratory 
200   Pleasant   Valley  V/ay 
West   Orange,   New   Jersey 

Westinghouse   Electric  Corporation 

Aviation   Gas   Turbine   Division 
iiTTN:      Engineering   Library 


P.  0.  Box  288 

Kansas  City,  Missouri 

FOREIGN  ORGANIZATIONS 


(1) 


Director  (1) 

National  Aeronautical  Establishment 
Ottawa,  Ontario 
Canada 


University  of  Toronto 
Institute  of  Aerophyslcs 
ATTN:   Library 
Toronto  St   Canada 


(1) 
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Trailing  Center  for  (1) 

Experimental  aerodynamics 
aTTN:      Library 

Rhode-Saint-Genese    (Belgique) 
72,   Chaussee   de   Waterloo 
Belgixim 

Chairman  (1) 

Defence   Research  Board 

ATTN:      DSIS 

Ottawa,    Ontario 

Canada 
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